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FOREWORD 

This  study  was  conducted  under  the  auspices  of  the  Defense  Mapping  Agency  in  order 
to  develop  formal  techniques  that  have  important  application  to  the  solution  of  propellant 
management  problems  related  to  the  selection  of  low-altitude  earth  satellite  flight  profiles. 
These  techniques  should  help  fill  the  analytical  vacuum  currently  associated  with  this  aspect 
of  mission  planning.  This  report  was  reviewed  and  approved  by  Dr.  R.  J.  Anderle  and  Mr. 
R.  W.  Hill. 


Released  by  ; 

T.  A.  CLARE.  Head 
Strategic  Systems  Department 
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INTRODUCTION 


Several  of  the  objectives  to  be  satisfied  by  earth  satellite  mission  planning  analysis 
are  to  predict  propellant  lifetimes  and  to  understand  the  various  propellant  conservation 
trade-offs  that  can  be  made  available  by  adopting  certain  flight  profiles.  This  aspect  of  mission 
planning  is  often  of  special  importance  for  low-altitude  missions,  where  considerable  amounts 
of  propellant  are  used  for  altitude  maintenance  thrusting  to  compensate  for  the  effects  of 
atmospheric  drag  deceleration. 

Two  methods  are  generally  employed  to  perform  propellant  consumption  studies; 

1.  The  application  of  empirically  derived  propellant  consumption  rates 

2.  The  utilization  of  analytic  approaches  based  upon  geopotential,  drag,  and  discreet 
thrusting  models  of  varying  complexity 

Both  of  these  methods  are  quite  effective.  However,  the  first  is  often  limited  by  the  lack  of 
data  describing  the  effects  of  varying  satellite  and  orbital  geometries,  as  well  as  by  solar 
and  atmospheric  conditions.  Tlie  second  approach  often  tends  to  be  inefficient  and  tedious  to 
use. 


Tliis  paper  presents  an  alternative  (and  perhaps  more  flexible)  analytic  development 
for  predicting  the  rate  at  which  propellant  is  consumed  for  in-track  drag  compensation  thrust¬ 
ing  for  low-altitude  earth  satellites.  The  associated  results  can  be  effectively  applied  to  pro¬ 
pellant  longevity  analyses,  especially  to  trade-off  studies.  Tlie  following  sections  discuss  this 
development  in  detail,  starting  with  derivations  of  the  propellant  longevity  and  mass  decre¬ 
ment  equations.  Tliese  are  then  applied  to  spherically  symmetric  atmospheres  that,  under 
certain  assumptions,  can  provide  analytically  tractable  results.  Such  results  are  useful  for 
comparison  with  those  obtained  from  more  complex  cases  under  special  limiting  conditions. 
Tlie  more  realistic  and  complex  case  of  propellant  consumption  in  an  oblate  diurnal  atmos¬ 
phere  is  also  considered.  An  analytic  expression  for  the  associated  mass  decrement  equation  is 
derived  and  used  in  the  following  section  to  provide  the  reader  with  several  illustrative  nu¬ 
merical  examples. 
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THE  PROPELLANT  LONGEVITY  EQUATION 


Consider  an  artil’ieial  earth  satellite  that  has  been  inserted  into  an  initial  nominal 
orbit  at  time  t^  with  a  total  weight  of  propellant  AW^  available  for  orbit  initialization  and 
maintenance  during  its  mission  life.  If  the  satellite  were  to  operate  in  N  different  nominal 
orbits  during  its  mission  life,  each  requiring  the  weight  of  propellant  to  achieve  the 

i'^  nominal  orbit  and  using  propellant  at  the  rate  (R.  during  the  interval  j  while  in  the 

i'*'  nominal  orbit,  then  the  following  expression  may  be  written; 


Note  that  AW^  ^  may  result  from  orbital  transfer  thrusting,  post-launch  orbit  initialization, 
and  de-orbit  thrusting. 

Assume  now  that  the  propellant  consumption  rates  dC  can  be  separated  into  a  radial 
thrusting  rate  ,  an  in-track  thrusting  rate  dlf.  a  cross-track  thrusting  rate  (R^",  and  a  miscel¬ 
laneous  operational  maintenance  consumption  rate  then, 

dC  =  dl*"  +  d?‘  dl*^'  +  d?”^’  (2) 

I  I  i  1  I  ^  ' 


and  Equation  ( 1  )  becomes 


I 


-I-  d?^' 
1 


AW 


0  A 


(3) 


The  prime  on  the  second  summation  in  the  last  equation  denotes  that  no  propellant  is  used 
if  natural  drag  decay  is  used  to  establish  shorter  period  orbits  and  that  the  summation  should 
not  include  these  natural  decay  intervals.  Tlie  time  bounds  for  such  intervals  are  determined 
the  .Jrae  decay  rates  for  those  intervals. 


Hu.  report  is  concerned  only  with  the  rate  at  which  [iropellant  is  spent  for  in-track 
;i,i  thus,  tile  torm  of  l,(|uation  (3)  may  be  simplified  by  defining 
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so  that 


It  is  simpler  to  express  Equation  (5)  in  terms  of  the  propellant  mass  decrement  Am  (i.e., 
the  in-track  propellant  mass  expenditure  for  one  orbital  revolution).  As  will  be  shown  in  the 
following  sections,  this  quantity  is  easily  derivable  for  certain  assumed  conditions.  This  simpli¬ 
fication  is  introduced  through  the  following  substitution; 


t 

/ 


(5^!  dt  -»•  g  T'  Am..  An.. 

1  “  ^  ij  I) 

j=  I 


(6) 


i-  1 


where  the  summation  is  taken  over  M  groups  of  orbital  revolutions  in  the  i"’  nominal  orbit 
An.j  for  which  the  mass  decrement  Am.,  applies  and  g  is  the  gravitational  acceleration. 
Equation  (5)  then  assumes  the  completely  discreet  form 


iW,  -  *'£  bL  Am„  An,  -  AW„^^ 

i=  1  '  j=l 


I  = 

I 


(7) 


Ellis  equation  shall  be  referred  to  in  the  following  sections  as  the  propellant  longevity  equation 
(PEE),  since  it  can  be  used  to  compute  the  propellant  life  X  defined  by 


N  M 

X  -  V  V  An  .  T..  +  AX„ 

^  ^  IJ  u  n 

1-1  j=i 


(8) 


where  r.  is  a  nominal  orbit  period  for  An  .  orbital  revolutions  and  AX,,  is  the  total  time 

ij  •  IJ  D 

spent  during  natural  decay  phases. 


THE  MASS  DECREMENT  EQUATION 


Assume  that  an  artificial  earth  satellite  in  a  nominal  orbit  with  semimajoj-  axis  a  and 
eccentricity  e  is  continuously  experiencing  atmospheric  drag  deceleration  and  is  sinuillaneously 
performing  in-track  microtlirusting  to  offset  the  effects  of  drag  decay  so  that  the  nominal 
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orbit  is  maintained.  Let  dE^  be  an  infinitesimal  orbital  energy  change  induced  by  drag  decay 
and  dE.j,  be  an  infinitensimal  energy  change  produced  by  the  microthrusting.  In  order 
that  the  nominal  orbit  be  maintained,  the  following  condition  must  be  satisfied: 

dE^  =  -  dEjj  (9) 

The  work  done  by  the  drag  force  (which  is  assumed  to  operate  only  in  the  direction 
opposite  the  satellite  velocity  vector  v  )  in  the  infinitesimal  time  element  dt  is 

dE|^  =  -  ^  C^Ap  lv|'  dt  (10) 

where  is  the  drag  coefficient,  A  is  the  cross  sectional  area  of  the  satellite  in  the 
direction  of  motion,  and  p  is  the  atmospheric  density  at  the  satellite’s  location.  The  6  factor 
accounts  for  the  orbital  orientation  with  respect  to  a  rotating  atmosphere  and  has  the  form 

where  r^  and  v^  are  the  radius  and  velocity  of  the  satellite  at  perigee,  respectively,  is 
the  earth’s  angular  rotation  rate,  i  is  the  orbital  inclination,  and  A  is  the  ratio  of  the  atmos¬ 
pheric  to  earth  angular  rotation  rates.  The  work  done  by  in-track  microthrusting  in  the  same 
infinitesimal  time  element  dt  is 

dE.^  =  T  |v  I  dt  (12) 

or.  alter  applying  the  “rocket  equation”, 

di:.^  =  -  g  l^p  iv  1  dm  (13) 

where  I  is  tile  thrust,  is  the  thruster  specific  impulse,  and  dm  is  the  infinitesimal  element 
‘>1  mass  of  the  proiiellant  e.xpended  during  the  thrust  interval  dt. 

Use  of  1  (piations  (10)  and  (13)  in  Equation  (9)  allows  the  in-track  propellant  ex- 

pmuhUire  rate  equation  for  drag  deceleration  offset  to  be  written  in  terms  of  the  satellite’s 

orbital  and  drae  eharaeteristics,  as  well  as  atmospheric  and  thruster  properties: 


4 
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dm 

dt 


A  S 


2g  I 


|V 


sp 


( i4) 


The  in-track  mass  expenditure  during  one  orbital  revolution  (i.e.,  the  mass  decrement)  is 
obtained  by  integrating  the  last  equation  over  one  orbital  period. 


Am 


2g  1 


r 


!V  l^p  dt 


(15) 


It  is  convenient  to  change  the  independent  variable  from  t  to  1:,  the  eccentric  anomaly. 
This  IS  done  using  the  following  relationships: 


Iv 


2 


+  ecosE 
-  e  cos  F 


and 


dt 


(1 


ecosE)dE 


(16) 


(17) 


where  p  is  the  earth’s  gravitational  constant.  Using  Equations  (16)  and  (17)  in  Equation  (15) 
and  making  the  upper  integration  bound  consistent  with  the  new  independent  variable  gives 

r„  A  6 

Am  =  -  — -  (pa) 

2e  I 

-  ,p 

In  (orinimj  tliis  expression,  it  is  assumed  that  ( .  A.  and  1  are  constant  over  the  re\o- 
lution.  rills  eciiiation  will  he  referrerl  to  hereafter  as  tltc  mass  decrement  equation  (.\1))l  I 


/ 


(1  +  eeosE)p  dE 


(18) 


PKOIM  LL\NT  CONSUMPTION  IN  SPHERICALLY  SYMMLTKU  AL 

ATMOSPHLRLS 


lee  li  i||i  1  iv  in.',  sul'sei  lions  .nx  devoted  to  the  consider.iiion  oi  liu  oiscs  wlrre  Ihe 
I  le' spi  :  h.  (!e;i-it\  Is  .issniiied  to  i'c  sphencaliv  symmetriCvd  riiesc  eases  ,iu  ol  mtie.t 
)  IS  n  Iheir  iiialvlie  tractahility  .Allhoueli  then  utilitv  ma\  'le  M  in/.vlni  li;e,iie.l  i  i  ', 

,n,  I m  vid  ml  oe-t'ii'j  msiglits  into  the  more  ecner.il  prold  -m 
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CONSTANT  DENSITY 


Tliis  subsection  is  concerned  with  the  case  where  the  atmospheric  density  is  assumed 
to  be  spherically  symmetrical  and  constant  with  altitude  above  the  surface  of  the  earth.  If 


P  ""  P,. 


where  p  is  a  constant,  then  Equation  (18)  becomes 


A  6 

Am  =  -  — -  (pa)"^  p^ 

^sp 


(1  +  ecosE)dE 


rrC^  A  5 


(pa)’'''  p„ 


This  result  may  be  readily  applied  to  the  PLE  to  give 


n'  (  ttC^ 

+E  I - T 


A  6  M  j 

-  PoT,  An..  -  1=  0  (22) 

p  j=i  \ 


Further  simplification  may  be  introduced  by  restricting  this  application  to  orbits  having  the 
same  nominal  orbital  period  (i.e.,  the  same  a),  but  different  eccentricities.  Since 


£  ^",1 

i-  1 


i:  n, 


V'.!  r'  n'  IS  tlie  Iol:il  number  of  orbital  revolutions  spent  in  the  N  nominal  orliits  (exchitiing 
'h  ~  i  cuirrmc  during  natural  tlecav  phases),  then  it  is  found  that 
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-  E  I 


TT  C,^  A  6  (Iiay^ 


I'hj  associated  propellant  lit'e  £  is  obtained  by  iisine  lu|Liatioiis  (23i  thrcrtieh  (25)  and  the 
relation 


^('  that 


-  -  E  a 


c;,,  A  5  p  p^, 


+  A£„ 


Tilts  result  substantiates  what  one  expects  intuitively  about  propellant  longevity,  i.e.. 
’fopellant  life  is  eninincetl  when 

1,  Less  propelhint  is  used  for  operational  maintenance  and  configuration  change 
i.niht  ailjusts 

2  The  nominal  semimajor  axis  is  increased 
fhe  thruster  specific  impulse  is  increased 
■i.  Hie  tlrag  ctK-fficient  and  satellite  cross  sectional  area  are  decreased 
The  atmospheric  density  is  decreased 


XI’OM  \  ri ALLY  DlC  RLASfNti  DENSITY 


In  tin  •'..inm,  an  ;it mos|ihcric  density  moilel  is  used  that  assumes  that  the  dcnsit\ 

I  I  ',  in’-'in  til,’  dotance  r  trom  the  center  of  tlie  earth  and  \anes  e\(ionentiall\ 
L  1^  d.'ioO',  model  ini'  the  torm 


III 


( 2  S  I 
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r,,  Is  til.'  I'.MCC  ilistaiK\.'  from  tlio  oartli'^  center,  p,,  is  the  i.lc:isil\  ai  ili. 

*  O 

aiiu 

! 

i 

J  —  ■ 

11 

‘-e.iv  ii  the  deiisitc  scale  height  and  is  assimied  constant.  .Siihslitntme  the  niai.i'ii- 


•■.(  1  -  eeosl  ) 


r,.  -  all  --  e) 


I  ■  i  Li  i tm;;  (  fS  i  ei'^  es 


e'\p  !  ■■  dae  +  pae  eosh.  | 


le:  tias  eepiessiar,  m  the  .\!l)h  gives 
(  ,  .\  6 


Air, 


2rr 

/ 


(p..(  '  pj,  e\p  (-  dsiel  1  (1  +  ecosl:)e\p  (due  cosl  Idl 


r(  ,  ,\  A 


Ip. 


Pi  e\p|-,)ae|  '  I  dpaei  +-  e  I^tpariJ 


.1  fii 


:;l  IP 


i'-^  I  «.  *  <  I  P.M  i  1' ' 
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t  ill  ai'.uii  i34)  can  lie  subslitiited  into  the  PLH  to  ^ive 


AW 


AW 


N  I  M 


N  '  1  A6 

-t-  ^  -  -  (jLtaA'-'-p  e\pt-j3.a.iv  1  )  I  ()3.a.e  ) 

A—/  i  f  I  '  ^1  11  I  1  1 


|  -  1 


I 

sp 


M 


J-  I 


AW 


0  A  . 


( 


0 


(.O' 


tl  ’.''-  siinplif'ieation  of  this  e\|iression  is  liitTicult  except  for  the  ease  where  N  =  1 .  Iheii. 
le  i 'Uiation  (24).  one  finds  that 


(AW,  -  AW^,^,  -  AW^J  1^,,  exp  [^ae] 
~  C,j  A5(, ua)'‘  p,,  ll  _(fjae)  +  e  l,(p'ae)l 


that  the  IM.f  heeontes 


:i^w, 


C  = 


AW  -  AW 

N  1  .\i  ^  o.\ 


a  exp(j3ae: 


C'l,  A  6pp,,  (l^^(j3ae)  +  e  I,(i3ae)l 


AX.. 


(38) 


expa■es^lon  also  suhstaiitiates  the  five  points  concerning  propellant  life  enhancetnenl  that 
iliscussed  earlier.  It  shonkl  be  noted  that  ficpiation  (38)  reduces  to  the  result  o; 
iiiini  (.'")  for  \  =  I  when  the  eccentricity  is  zero,  since 

I  I  i  I  I  .  i  (.Oil 


i’ROf’I  LLXNT  COiN .SUMPTION  IN  AN  OBLATt  ATMOSPHKRi 
WIfU  DAY-iO-NIOHT  IM  NSITY  VARlAllON 

eiii';i,  a)i  an.ilylie  exjiression  is  developed  lor  the  Ml)!  tor  a  low  altit:,. 

'  i'  '..tiri'’  :n  ,):i  orbit  with  small  eeecntrieity  (0.01  •_;  e  C  O.liO  1  he  lonu  oi 
h  ;re,t,  model  is  tlial  of  an  olslate  atmosphere  witli  d,i\ -I o-iiieli t  uuisi! 


NSWC  TR  83-243 


An  analytic  form  for  the  atmosplieric  density  can  be  obtained  by  conihininy’  the 
oblate  atmosphere  I'onn  used  by  Cook,  King-Hele,  and  Walker’  with  that  of  an  atmo^idiere 
with  diurnal  variation  discussed  isv  Cook  and  Kine-Hele^.  The  resulting  torm  lor  the  densit\ 
is  given  b\ 


p  =  (1  +  Fcos0)  exp  1-  p'ir  -  a)|  (40) 

where  p.  is  a  reference  atmospheric  density  of  the  form 


1 

p  -  -  ip  +  p  .  ) 

o  1  max  ni  in 


P  -  P  ■ 

ni  a  X  m  I  n 


F  =  - 

P 

^  Ill  a  X 

+  p 

min 

(42) 

*7 

It 

(43) 

^  1 

-  e  sin^i  sin^Li  ( 

(44) 

(7  —  r  *1 

1 

-  €  sin^i  sin^co  \ 

and 


COS0 


cos  T  -  e 
1  -  e  cos  1: 


+  B 


( 1  -  e^ sin  T 
1  -  e  cos  K 


(4.^) 


In  the  last  four  equations,  p  and  p  are  the  maximum  day  time  and  minimum  nieht 
time  den.xities,  respectively;  H  is  the  density  scale  heiglit;  i  and  co  are  the  orbital  inclination 
and  argument  of  perigee,  respectively;  u  (=  co  +  0.  where  0  is  the  true  anomaly)  and  T  are 
tile  true  argument  of  latitude  and  eccentric  a  omaly.  respectiveh ;  e  is  the  earth's  ellipticitc'; 

IS  the  perigeal  ratlins;  and 


A  -  xin  sin  i  sin 


CO  +  cos  ;  cos(12  -  )cos  CJ 


ctis  i  sin  (i2 


)  sin  CO 


(4i) ) 


and 


H  sin  6 sin  i  cos  CO  -  cos  h  ■  cos 


sin  CO  +  cos  I  sm(i2  i  cos 


CO 


(40 


'  (,  I  I,  Kiti  >:  I  i  i' Ic  ,  iin  a  Wiill.iT.  /V>i  I  ('.//nv'  '>(  I’/ic  Woi  u/ ,SVi(  7(71  1 .  2(>4 ,  p  r  s  s  1  :  I .  I 'ii,  i , 

^  ,  J  f  . .  I  V  ,  in  1  I  >-  O  K  up:  I  in  Ic.  /*rf  ui' t'J  nms  *>]  ihr  AC  M"  7,7  .S',  .1,  2.S9.  pp.  .1.1  f,7.  I  ‘0.5  . 


NSWC  TR  83-243 


where  and  are  the  right  ascension  and  declination  of  the  atmosplieric  diurnal  bulge, 
respectively,  and  S2  is  the  right  ascension  of  the  ascending  node  of  the  satellite  orbit. 

liquation  (44)  may  be  expanded  to  first-order  in  e  to  give 

CT  =  r^  t  1+1/26  sin^i  (cos2u  -  cos2cl))]  (48) 

Similarly.  Equation  (45)  may  be  expanded  to  first-order  in  e  to  give 

cos  <p  =  A(cosE  -  e  +  ecos^  E)  +  B(sinE  +  csinEcosE)  (49) 

Substituting  Equations  (48)  and  (49)  into  Equation  (40)  and  using  the  relation 

r  =  a(l  -  ecosE)  (50) 

allows  the  following  first-order  expression  to  be  written  for  the  atmospheric  density; 

fj  =  p  [1  +  FA(cosE  -  e  +  ecos^E)  +  FB(sin  E  +  e  sinE  cosE)] 

(51) 

exp  j3ae  (1  -  cosE)  +  c  cos2u  -  c  cos  2  cu  | 

vv  here 


c  =  -  6  i3  r^,  sin^i 


(521 


As  discussed  in  Retcrence  1,  c  may  be  treated  as  a  small  parameter  of  the  same  order  ot 
ni.i'jnituile  as  the  eccentricity.  Ihus.  in  lupiation  (51),  the  following  expansion  may  be  used 


exp  '  c  cos 


+  c  cos  2 11  + 
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Substituting  liquations  (51)  and  (53)  into  Equation  (18)  and  using  the  relations 


eos  E  - 


-  e^)'''-  sin  E 


to  cliniinate  0  to  first  order  in  e  gives 


Am  =  -  — - exp  {-  ^ae  -ccos2a))  I  ill  FA(cos  E  -  e  +  ecos^E) 

sp 


+  FB(sin  E  +  e  siiiEcosE)]  •[  1  +  ccos2(cj  +  E)  -  2ec  sin2(co  +  E)  sin  E 


+  -  eos4(oj  +  E)  -  ec^  sin  4(c«;  +  E)  sin  E) 

4  4 


i  I  +  e  eos  E.  I  exp  [)3ae  cos  E)  |  dE 


When  the  integrand  ot  this  equation  is  multiplied,  the  result  contains  trigonometric  terms 
that  are  expressible  as  fiinctions  of  cos(nE),  n  =  0,  1,  2,  6,  Tliis  allows  Equation  (56) 

to  lie  written  in  terms  of  the  integral  representation  of  the  Bessel  function  of  the  first  kind 
and  imagmarv  argument  defined  by 


(e!ae)  =  —  I  cos(nI-j  exp((3ae  cos  E)dE 

'  2vr  f 


i  '.ultme  MDl  lor  an  oblate  diurnal  atmosphere  is 
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Am  = 


TT  Cp  A  6 


g  I 


sp 


p  exp  I  -  Pae  -  ccos2  cj  | 


(-t) 


1  +-j  (I^  + 


FA 


t) 


1  +  —  )  (I,  +  el^) 


+  j  {  [21^  -  e(I^  -  3I3)]  +  FA  [I,  +  I3  + 


:elj  }  cos2cj  -  ^FB  j  (I,  -  I3)  +  2e  (I^  -  I^)(  sin  2  w  + 


(58) 


{[21^  -  e(3l3  -  5I3)]  +  FA  [(I3  +  Ij)  +  e  (31^  -  I^)]  }  cos  4  w  + 


.‘2 


—  FB  {  [I5  -  I3  ]  +  €[(3  -  41^  +316^  sin  4  co 


where  the  Bessel  function  argument  0ae  has  been  suppressed  for  the  sake  of  brevity. 

Tliis  MDF  could  be  introduced  into  the  PLE  at  this  point.  However,  due  to  its  com- 
j  plexity  and  dependence  upon  solar  position  and  the  additional  orbital  parameters  n  and  u>. 

little,  if  any,  simplification  could  be  introduced  into  the  resulting  analytic  expression.  It 
should  be  noted,  nonetheless,  that  Equation  (58)  reduces  to  the  form  of  Equation  (34) 
when  oblateness  and  diurnal  effects  are  neglected;  i.e.,  when 

I  A  0 

B  ^  0 


and 


J 


c  -  0 


NUMERICAL  EXAMPLES 

Several  numerical  examples  that  illustrate  the  types  of  analyses  to  which  the  MDF  i'v 
an  oblate  diurnal  atmosphere  may  be  applied  are  presented  in  this  section.  These  examples 
were  created  using  analytically  modeled  averaged  variational  equations^  to  represent  the  effects 

I.  J.  I  I. ill  an>l  K.  L.  AlfVtrd.  "Se m ia n a ly I ic  Theory  for  a  Close-F^arth  Satellite,"  /1//1/1  Journal  of  (inidancc  atui 
intitrol.  Void.  Ju  Iv  Aiinust  I  h  RO.  p  p  .  .104  3  1  I . 
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ot  gcopotoiitial  perturbations  on  the  satellite  motion  through  J4.  Theoretical  expressions 
for  drag  decay  rates  for  a  low-altitude  satellite  orbiting  in  an  oblate  diurnal  atmosphere^  were 
used  to  model  the  natural  drag  decay  phase  occurring  between  the  time  of  propellant  de¬ 
pletion  and  reentry.  Tlie  MDE  of  Equation  (58)  was  used  to  compute  the  propellant  re¬ 
quirements  for  orbit  sustenance.  Changes  in  orbit  parameters  occurring  during  orbit  adjusts 
were  computed  using  results  obtained  from  the  Lagrange  planetary  equations  when  impulsive 
velocity  clianges  are  assumed.^  The  quantities  of  propellant  expended  during  an  orbit  adjust 
were  obtained  from  application  of  the  “rocket  equation.” 

The  Jacchia  1960  model  atmosphere  was  used  for  all  density  computations  required 
for  both  the  MDE  and  the  drag  decay  rates:  i.e.,  computation  of  the  density  at  the 

diurnal  bulge  location  and  the  density  diametrically  opposite  the  bulge  at  the  satellite’s 

osculating  perigee  altitude.  Tliis  model  atmosphere  describes  the  density  variation  with  altitude 
of  an  oblate  diurnal  atmosphere  and  accounts  for  the  effects  of  density  variation  due  to  solar 
activity  via  a  dependance  upon  the  solar  flux  F,q7.  It  is  not  believed  to  be  extremely 
representative  of  the  atmospheric  density,  but  is  computationally  very  efficient.  Santora's  method 
for  density  scale  height  selection**  was  employed  for  both  drag  decay  rate  and  MDE  evaluation. 

The  following  data  were  used  to  initialize  the  computations  for  each  of  the  examples 

below : 


a  =  6756.205  km 
e  =  0.009656113 
r  =  94.99996° 
w  =  193.3874° 

{1  =  95.00015°  (59) 

M  =  156.61  15° 

i,_  =  44619.987  Modified  Julian  Days 
1,,,  =  230  sec 
AW,  =  150  kg 

AVI  R  \(,|  IN  TRACK  CONSUMPTION  RATE  PREDICTION 

f  Dinputations  were  iierformed  using  the  conditions  of  Eciuation  (59)  for  F,  g  7  values 
I'  lwvcii  100  and  300  Dux  units.  Three  C„A  values  of  2.  x  10“'*  km^.  1.  x  10'*  km^. 
and  \  10  ^  km^  were  assumed.  The  average  in-track  consumption  rate  (O'  was  obtained 


I  \  s*tit'<ri  ’  ^.tt  Miff  Dr.i.k:  I’frturhal  ions  in  an  OblaCe  Diurnal  Afmosphere,*' /I //I /I  Journal.  Vol.  I  3,  Scptem her  1975,  pp.  121 
! :  1 1- 

\  h  1  /t>Kf,;tiii.'nrnus  Orhit  l^aramcicr  f '/w/jifcv  f*roiJuced  hv  /mpulsiicThrustinfiwith  Application  to  Orhit  Adjust  Dcsii;n  for 
'  x'o,;// /  M  ( /ffr/c/f  i  NSWC  /DL  I  K  83-31  (Dahl^jrt'n.  Va..  I  chuniry  1983).  pp  3  4. 
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lor  each  F,,,,  and  Cj^A  combination  by  dividing  the  150  kg  of  propellant  consumed  by  the 
time  required  for  its  consumption  (i.e.,  the  length  of  time  the  orbit  was  sustained).  The 
results  are  presented  in  Figure  1.  It  should  be  noted  that  those  rates  are  representative  only 
of  those  reciuired  for  sustenance  of  the  orbit  described  by  the  elements  of  Equation  (59). 
since  prereentry  drag  decay  phases  are  not  included. 


FIGURE  1.  AVERAGE  IN-TRACK  PROPELLANT  CONSUMPTION 
RATES  AS  A  FUNCTION  OF  SOLAR  FLUX  CON¬ 
DITIONS  AND  SATELLITE  DRAG  CHARACTERISTICS 


SAMPLE  FLIGHT  PROFILE  TRADE-OFF  ANALYSIS 


Consider  tlie  following  scenario:  “A  satellite  mission  is  ordinarily  constrained  to  operate 
m  an  orbit  described  by  the  parameters  of  Equation  (59)  until  its  propellant  is  depleted. 
It  whicli  time  it  deorbits  by  natural  drag  decay.  After  orbital  insertion,  it  becomes  apparent 
that  It  will  likely  be  necessary  to  take  measures  at  sometime  during  the  mission  to  extend 
the  mission  length  for  as  long  as  possible.”  The  MDE  can  be  applied  to  the  situation  de¬ 
picted  by  this  scenario  to  provide  estimates  of  the  propellant  and  orbital  lifetime  trade-offs 


th.il  exist. 
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kvsLlll^  aix'  |'i\M.'iUod  lor  tlic  spL'cial  case  where  tlie  semiinajor  axis  is  increased  30  kin 
In  .111  i!i-track  apoeee  thrust  at  some  time  iliirint:  tlie  mission  and  tlie  satellite  deorbits  by 

iiaiui.d  diae  deeav.  A  representative  tliplit  profile  is  presented  m  Figure  2.  Ti  desired  data 

neia  eeiiei  ited  Using  1  ,  „  ,  =  100  and  (  A  =  I.  \  1 0" km^  and  are  shown  in  Figure  3. 
llui.  ihe  change  m  the  propellant  longevitv  and  the  increase  in  total  mission  life  are  plotted 
,1-Muni  ilic  icvolution  number  during  which  the  orbit  aiijust  was  perlornied.  F'or  example. 

:i  !:ic  aiiiusl  was  pertomied  on  rev  420.  there  would  be  no  change  in  the  propellant  longevity 
I'.ci  ti,al  il  no  adjust  was  pertormevl.  I'he  mission  life  would  also  be  extended  by  190  days 

dn.'  Ill  the  decrease  m  the  drag  decay  rate  obtained  by  operating  at  a  higher  altitude. 


FKiURL  2  RLPRESLNTATIVL  FLIGHT  PROFILE 
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FIGURE  3.  IMPACT  OF  ORBIT  MODIFICATION  UPON 
PROPELLANT  AND  MISSION  LIFETIMES 


SUMMARY 


Tlie  propellant  longevity  equation  has  been  derived  in  both  its  continuous  and  discreet 
forms  and  expressions  for  the  mass  decrement  equation  in  spherically  symmetrical  and  oblate 
diurnal  atmospheres  have  been  developed.  Tlie  mass  deerement  equation  tor  an  oblate  diurnal 
atmosphere  has  been  applied  to  the  discreet  I'orm  of  llie  propellant  longevity  equ.ilion  to 
provide  several  numerical  examples  that  illustrate  their  applicability  to  the  solution  ot  certain 
mission  planning  problems. 
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